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Abstract 

In this paper we construct nonlinear partial differential equations in more than 3 in- 
dependent variables, possessing a manifold of analytic solutions with high, but not full, 
dimensionality. For this reason we call them "partially integrable" . Such a construction is 
achieved using a suitable modification of the classical dressing scheme, consisting in assum- 
ing that the kernel of the basic integral operator of the dressing formalism be nontrivial. 
This new hypothesis leads to the construction of: 1) a linear system of compatible spectral 
problems for the solution U(X; x) of the integral equation in 3 independent variables each 
(while the usual dressing method generates spectral problems in 1 or 2 dimensions); 2) 
a system of nonlinear partial differential equations in n dimensions (n > 3), possessing a 
manifold of analytic solutions of dimension (n — 2), which includes one largely arbitrary 
relation among the fields. These nonlinear equations can also contain an arbitrary forcing. 



1 Introduction 

Since the discovery of the integrability of the Korteweg-de Vries equation [1], much effort 
has been devoted to the study of direct techniques to construct and solve nonlinear Partial 
Differential Equations (PDEs). One of the most powerful of such techniques is the dressing 
method, originally introduced in [2] and subsequently generalized in [3]- [6] (see also [7, 8]), 
which is based on the existence of a linear analyticity problem, i.e. a Riemann-Hilbert or a 8 
problem in some spectral variable A for some matrix function U(X; x), depending parametrically 
on the space-time variables X — (x X ). (The d problem was introduced, in the context of 
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integrable systems, in [9, 10].) Such an analyticity problem is characterized by a linear integral 
equation, whose unique solvability allows one to construct and solve an overdetermined system 
of compatible linear spectral problems for U(X;x), and, consequently, a nonlinear system of 
PDEs in the independent variables x, for the coefficients of such a linear system. 

The manifold of the analytic solutions of the nonlinear PDEs constructed by the dressing 
method is parameterized in terms of a proper number of arbitrary spectral functions, appearing 
in the linear integral equation, which depend onn-1 variables. Therefore the solution space is 
full and we say that the nonlinear PDE is completely integrable. For instance, the solution space 
of 1 + 1 dimensional scalar systems like the Korteweg-de Vries and the nonlinear Schrodinger 
equations has dimension 1, being parameterized by an arbitrary function of 1 variable, while 
the solution space of the 2 + 1 dimensional generalizations of them, the Kadomtsev-Petviashvili 
and the Davey-Stewartson equations, has dimensionality 2, being parameterized by an arbitrary 
function of 2 variables. 

Motivated by the above considerations, in this paper we say that the dimensionality of the 
space of analytic solutions of a system of PDEs is k, if the analytic solutions are parameterized 
by a "sufficient number" of arbitrary functions of k independent variables. For instance, if the 
system of PDEs contains K equations involving first order "time" -derivatives of K functions, 
then the "sufficient number" is K. If k = n — 1, then the systems is completely integrable. 

One of the most outstanding open problems in the theory of integrable systems is the 
construction of nonlinear PDEs in multi-dimensions (see, for instance, [12, 13]), i.e., in more 
than 3 dimensions, which could be integrated using suitable extensions of the above dressing 
procedure. Apart from few exceptional instances, among which one counts the self-dual Yang- 
Mills equations [11] and the Plebanski heavenly equation [14] (see [15, 16, 17, 18] and [19, 20] 
for their integration schemes), no significant examples are known in the literature [21]. 

The purpose of this paper is the construction of PDEs in more than 3 independent variables, 
possessing a manifold of analytic solutions with high, but not full, dimensionality. For this 
reason we call such PDEs "partially integrable" . This construction is achieved using a suitable 
modification of the classical dressing scheme, consisting in assuming that the kernel of the basic 
integral operator of the dressing formalism be nontrivial. As we shall see, this new hypothesis 
leads to the construction of: 

1. a linear system of compatible spectral problems in 3 independent variables each, for the 
eigenfunction U(X;x), where A is a vector spectral parameter (while the usual dressing 
method generates spectral problems in 1 or 2 dimensions with scalar spectral parameter); 

2. partially integrable nonlinear PDEs in n dimensions, possessing a solution space of di- 
mension (n — 2). 

A prototype example is given by the following 4 dimensional system of two matrix equations 
B 2 (q 1 ,q 1 ,q 2 )B 2 ~ 1 (q 1 ,q 2 ,q 3 ) = B 3 (qi, qi, q 2 )B^ 1 (q 1 , q 2 , q 3 ) = B 4 (qi, qi, q 2 )Bl 1 (qi, <?2, qs) (1) 
for the three square matrix fields qi(x), q 2 (x), qz{x), supplemented by the "largely arbitrary" 



relation 



F(qi,q2,qs) = 



(2) 



among them, where the matrix blocks Bj are defined as: 



Bj(qi, q 2 , q 3 ) = fei - fei-Bj - 55 Bj] - [Bj, q 3 ], j = 2, 3, 4, 



(3) 
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Bj, j = 2,3,4 are constant diagonal matrices different from the identity, and [•, ■] is the usual 
commutator of matrices. In the simplest case, the largely arbitrary relation (2) can be chosen 
to be an equation defining one of the fields, say q 3 , to be any given function j(x) (in general, 
a generalized function), interpretable as an "external arbitrary forcing": 

F: q 3 (x)= 1 (x). (4) 

As we shall see in the following, equations (1-4) possess a manifold of analytic solutions of 
dimension 2. 

The above closed system of equations (1-2) share with the other models constructed in this 
paper the following properties. 

1. The existence of a nontrivial kernel of the basic integral equation implies that the solutions 
constructed by the dressing depend on an arbitrary function f(x) of the coordinates; this 
fact has the following important implications. 

2. The nonlinear system of PDEs constructed by the dressing scheme possesses a distin- 
guished block structure (see (1)) and is under-determined. 

3. To close the system and fix its indeterminacy (or, equivalently, to fix f(x)), one has 
to introduce an "external and largely arbitrary" relation among the fields (see (2)). If, 
for instance, such a relation (algebraic or differential) is linear, then the construction of 
explicit solutions via the dressing algorithm remains linear as well. The simplest example 
of linear relation is obtained imposing that one of the fields be a given function of the 
coordinates, like in (4), interpretable as an external forcing. 

4. The system of PDEs depends on two types of matrix fields, those obtained "saturating the 
vector parameter A" of the solution U (A; x) of the linear integral equation by ingredients 
of the classical dressing method, whose dimensionality is constrained, and those obtained 
saturating A by a novel dressing function G(A; x), whose dimensionality is not constrained. 
That's why the dimensionality of the solution space, (n — 2), can be arbitrarily large. 

5. While integrable PDEs in low dimensions (2+1 or less) are the compatibility of overdeter- 
mined systems of linear spectral problems, such a feature seems to be lost for our higher 
dimensional examples. 

We remark that partially integrable equations of the type (1) are somehow connected to the 
N-wave type systems; indeed the two equations 

B 2 (qi, qi, q 2 ) = B 3 (qi, qi, q 2 ) = (5) 

are equivalent to the N-wave system in 2+1 dimensions for the field q 1 , obtained eliminating q 2 
from equations (5) (see also Sec. 2.1). The construction of partially integrable PDEs connected 
to other basic integrable systems, like those belonging to the Kadomtsev-Petviashvili hierarchy, 
or those associated with the Davey-Stewartson equation, will be the subject of a forthcoming 
paper. 

We also remark that a different generalization of the dressing procedure, allowing to con- 
struct a class of partially integrable PDEs which combine S and C integrability features, has 
been already proposed in [22] (a nonlinear system is S integrable, like the KdV equation, if 
it is solved via a linear integral equation; it is C integrable, like the Burgers equation, if it 
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is linearized by a simpler change of variables, like a contact transformation; see [23] for more 
details on these definitions). 

The paper is organized as follows. In Sec. 2 we derive, for the sake of comparison, the 
well-known iV-wave system in (2+l)-dimensions, using the classical dressing method, and we 
show that its solution space is full. Then in Sec. 3 we explore the implications of the existence 
of a nontrivial kernel of the basic integral operator of the dressing scheme, and we construct 
examples of nonlinear n-dimensional PDEs possessing a space of analytic solutions of dimension 
n — 2. We also show that these equations do not seem to be the commutativity condition of 
overdetermined systems of linear PDEs. In Sec. 4 we show how to construct an integral operator 
possessing a nontrivial kernel, and we use this result to characterize a large class of analytic 
solutions of the partially integrable PDEs of this paper. A natural extension of the algorithm 
presented in this paper is briefly mentioned in the final Sec. 5. 



2 Derivation of classical integrable systems using the 
dressing algorithm 



The starting point of the dressing constructions contained in this paper is the linear integral 
equation 



in the spectral variables A, fi, for the unknown matrix function U. The given matrix functions $ 
and \1/ are defined by some extra conditions, which fix their dependence on an additional vector 
parameter ), whose components are the independent variables of the associated 

nonlinear PDEs. Q is some largely arbitrary scalar measure in the /z-space. Apart from Q, all 
the functions appearing in this paper are Q x Q matrix functions. 

We remark that no a priori assumption is made in (6) on the dependence of ^ on A (this 
general starting point has been used, for instance, in [24] and in [22]), to keep the structure 
of \I> as much general as possible. Indeed, although in most of the cases such a dependence is 
described by a Cauchy kernel, an indication that equation (6) is the manifestation of Riemann- 
Hilbert and/or 8 analyticity problems, there are examples (see [24] and [22]) in which more 
general representations appear, indicating that the above analyticity problems could be a too 
restrictive starting points. 

Before developing, in Sec. 3, 4, the novel features of the dressing method, it is useful to 
summarize the essential steps of the classical dressing method used to construct and solve the 
classical 3-dimensional iV-wave system (18) (which is known to be an S'-integrable system), 
together with its solution space. Such a solution space is 2 dimensional (i.e., it is complete), 
being parameterized by an arbitrary function of 2 variables. 

2.1 ^-integrable PDEs: the iV-wave system. 

The basic assumption underlying all the known dressing procedures available in the literature, 
is that the operator ^ in (6) be uniquely invertible; i.e., that 



$(A;x) 



/ 



^(A, yu; x)U(fx; x)dQ(fj,) = (A; x), 



(6) 



dimker\l/ = 0. 



(7) 
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The ^-dependence is introduced by the matrix equations 

ty x i(\, n\ x) — <3>(A; x)BiC(fj,; x), i — 1, .., dim x, (8) 

showing that the x-derivatives of the kernel \1/ are degenerate matrix functions of the spectral 
parameters, another basic feature of all known dressing algorithms, where B i} % — 1, ..,n, are 
constant diagonal matrices, so at most Q of them may be independent. Due to the above 
degeneracy, the compatibility of equations (8) leads to separate equations for <3> and C: 

'I\ l>j 'I>, 1>, 0. i^j, (9) 

a/', r>,c, o, i^j, (io) 

and one equation is the adjoint of the other. Without loss of generality we assume B\ = I, 
where I is the identity matrix. 

Replacing, in equation (9), $ by as indicated in (6), and using (8), one obtains the 
following equation: 

VL tJ U = 0, (11) 

where 

LijU = U x iBj — U x jBi + UBivBj — UBjvBi, i,j = 1, .., dim x, i ^ j (12) 

and 

v{x) = J C(X;x)U{X;x)dQ(\). (13) 

Then the property (7) implies that 

LijU(X; x) = 0, i,j = 1, ..,dim x, i^j (14) 

or, explicitly: 

L 21 U = U X 2 - U x iB 2 - U[v, B 2 ] = 0, (15) 
L 31 U = U x s - U x iB 3 - U[v, B 3 ] = 0, 



having chosen j — 1, i — 2, 3. 

This is nothing but the well-known linear overdetermined system corresponding to the iV- 
wave equation in the three variables 

The associated complete system of nonlinear PDEs is simply obtained, in the dressing 
philosophy, upon "saturating the parameter A" in equations (15) by the integral operator 
JdQ(X)C(X;x)-. 

L 21 v - [B 2 , v 1 } = v x 2 - v x iB 2 - v[v, B 2 ] - [B 2 , v 1 ] = 0, (16) 
L 31 v - [B 3 , v 1 ] = v x s - v x iB 3 - v[v, B 3 ] - [B 3 , v 1 } = 0. 

It is written in terms of the square matrix fields v(x) and v 1 (x), where 

v\x) = J C x i(X;x)U(X;x)dQ(X). (17) 
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Eliminating v 1 from these two equations, we get the celebrated iV-wave system in 3 dimensions: 



[v x s, B 2 ] - [v x2 ,B 3 ] + B 2 v x iB 3 - B 3 v x iB 2 - [[v, B 2 ], [v, B 3 ]] = 0. (18) 

The same equation may be derived directly from the compatibility condition of the system (15). 

Similarly, considering the equations Lj-JJ = and L kl U = for any j ^ k ^ 1, one derives 
the hierarchy of n-wave equations 

[v x k,Bj] - [v xj , B k ] + BjV x iB k - B k v x iBj - [[v, Bj], [v, B k ]] = 0. (19) 

We remark that, in the above dressing construction, the linear integral operator \l/ in (6) 
acts from the left and, consequently, the partial differential operators in (14) act from the 
right, while, in the soliton literature, one usually makes the opposite choice. 

Our choice is motivated by the fact that, as we shall see in Sec. 3, 4, in more than 2 + 1 
dimensions, the role played by the linear integral equation (6) seems to be more fundamental 
than that played by the linear overdetermined system of PDEs. Indeed, while integrable PDEs 
in 2+1 dimensions (or less) are characterized as the compatibility condition of a linear overde- 
termined system of PDEs, such a basic property seems to be lost in multi-dimensions. Instead, 
as we shall see in the following sections, the linear integral equation (6) can generate nonlinear 
PDEs, together with their large manifold of analytic solutions, also in the multidimensional 
context. 

2.2 Solution space 

We now consider the manifold of particular solutions of equations (18,19). The solutions of 
eqs.(9) and (10) can be parameterized as follows: 

$(A;x) = J %(X,k)e kB - x dk, (20) 
C(fi;x) = J e qB - x C (fi,q)dq, (21) 

where 

n 

B ■x = ^B i x\ (22) 

i=l 

and where the spectral parameters A, /x, k, q are necessarily scalars. Thus equations (8) yield: 
*(A, Kx) = j $ (A, k)e^ B -*C (», q)^ + E(A, p), B 1 = I, (23) 

where the integration constant S(A, fi) is chosen here to be S(X — //). 

It is simple to see, from the linear limit, that the solution space of equation (18), generated 
by the dressing algorithm, is full. Indeed, in the linear limit: \I/(A, n) ~ S(X — ji) and U ~ <E>. 
Take C (A, q) = 8{X — q); then the solution v of the 3 - dimensional A^-wave system (18), which 
in the linear limit reads 

v{x)~ J C(X;x)<5>(X;x)dQ(X) = J e XB x %(X,k)e kB - x dkdQ(X), (24) 
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is parameterized by the arbitrary matrix function $o(A, k) of the two scalar spectral parameters 
A, k; then its solution space is 2 dimensional, and therefore it is complete. 

We end this section remarking that the Cauchy kernel appearing in (23), obtained here as a 
consequence of equations (20), (21) and (8), is a manifestation of the distinguished analyticity 
properties of the solution U(X;x) in the complex A plane, in agreement with the well-known 
derivations of the A^-wave equation (18) from Riemann-Hilbert and /or 3 problems [25, 26, 27]. 

3 Partially integrable PDEs in multi-dimensions. 

In this section we show how to construct partially integrable PDEs in n dimensions exhibiting 
a space of analytic solutions of dimension (n — 2). 

3.1 Generalization of the dressing algorithm 

In the Sec. 2.1 we have constructed, from the general hypothesis (6), (7) and (8) underlying the 
classical dressing algorithm, the integrable iV-wave system in 3 dimensions. The main obstacle 
to go to higher dimensions is clearly due to the fact that each linear problem L^U = 0, as a 
consequence of (7), is 2 dimensional. 

To increase the dimensionality of the linear problems, we then suppose that the kernel of 
the operator ^ is one dimensional: 

dimker^ = 1; (25) 

i.e., the solution of the homogeneous equation associated with eq.(6) is nontrivial: 

= ^H H(X;x) = U h (X;x)f(x), (26) 

where U h (X;x) is some nontrivial solution of the homogeneous equation ^>H = 0, f(x) is an 
arbitrary matrix function of x, and A, //, v are now vector spectral parameters whose dimension 
is specified in Sec. 4. Then the general solution of eq.(6) reads 

U(X-x) = U p (X;x) + U h (X;x)f(x), (27) 

where U P (X; x) is some particular solution of (6). 

As a consequence of the novel assumption (25), equation (11) implies the following equations 
for U: 

L l3 U{X-x) = {L nm U{X-x))A^ nm {x), i^j, n + m, (28) 

where A^ nm {x) are some matrix functions of x to be specified, reflecting the fact that two 
solutions of the homogeneous equation ^>H = differ by a matrix function of x, multiplied 
from right. 

Since the following cyclic permutation formula among three operators holds: 

£ (L tJ M)B k = 0, i^j^k, (29) 

cycl(ijk) 
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where M. is an arbitrary square matrix, it follows that only (n — 1) operators Ly are linearly 
independent. Therefore we take the operators {L jl: j — 2, . . . , n} as elements of the basis, and 
we consider the following subset of equations (28), involving only these elements: 

£j(\;x) = L jl U(X;x) - (L 2 iU(\;x))A j (x) — 0, j = 3, ...,n, (30) 
LnU^Uj-UsiBj-UfaBj], j = 2, . . . , dim x, (31) 

where A^rr) are some matrix functions to be defined. 

We have established that, if dimker\l/ = 1, then each linear equation (30) for the spectral 
function U (A; x) is 3 dimensional. 

The associated nonlinear equations, obtained "saturating the parameter A" in equations 
(30) by the integral operator j dQ,(X)C{X;x)-, read 

Ljiv - [B j ,v 1 ] - {L 21 v - [B 2 , v 1 ])^ = 0, 3 = 3, .., n; (32) 

they are given in terms of the fields v(x) and v 1 (x), defined respectively in (13) and (17), and 
of the matrices A^(x). More explicitly, one obtains: 

v xj - v x iBj - v[v, Bj\ - [Bj, v 1 ] - (v x 2 - v x iB 2 - v[v, B 2 ] - [B 2 , v 1 ]) A 3 =0, j = 3, .., n. 

(33) 

In order to express A^(x) in terms of U and close the system, we introduce an external 
dressing function G(X] x), and the associated new matrix fields 



w 00 {x) = I G(X;x)U(X;x)dn(X), w 3 °{x) = j G xJ {X; x)U{X; x)dQ{X), j > 0, (34) 
w lj (x) = J G x i xj (X;x)U(X;x)dQ(X), i,j>0, w ij (x) = w ji (x). (35) 



The equations for the fields w %i can be derived applying f dQ(X)G(X)- and f dQ(X)G x n(X)- 
to the linear equation (30), obtaining: 

L jlW n0 - w jn + w ln Bj = (L 21 w n0 - w 2n + w ln B 2 )A j , j = 3, .., dimx, n = 0, 1, .., dimx. (36) 

Some of these equations can be taken as definition of A> . But, to close the system, one 
needs to introduce the following additional structures. 

(a) Equations defining G(X;x). Since G is an outer dressing function, these equations can be 
a quite arbitrary system of compatible and solvable partial differential equations, either linear 
or nonlinear, involving derivatives of any order and dimension. This freedom plays a key role 
in allowing for a large solution space. 

(b) An additional relation between all the matrix fields, which may be taken in quite arbitrary 
form 

F(v, v 1 , w 00 , w i0 , w ij ) — 0, ij = 1,2, .... (37) 

This equation is needed to fix the arbitrary function f(x) of the variables x % appearing in the 
solution of the inhomogeneous equation (6) (see equations (26,27)). The relation (37) is largely 
arbitrary; the only requirement is that it should give rise to a solvable equation for f(x). Since 
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U depends linearly on the arbitrary function f(x), all the fields depend linearly on / as well. 
Indeed, using equations (13, 17, 34) and equation (27), one obtains 



v(x 
v 1 (x 



w 00 (x 



w i0 (x 



w lJ (x 



hl{x) + h\{x)f{x), h v (x = J C(X)U p (X)dtt(X), h\ = J C(X)U h {X)dQ(X), (38) 

h v \x) + hf(x)f(x), hf = Jc x i(\)u p (X)dn(\), hf = j c x i(x)u h (\)dn(\), 

hf°(x) + hf°(x)f(x), hf° = J G(X)U p (X)dQ(X), hf = J G(X)U h (X)dQ(X), 
hf(x) + hf(x)f(x), hf = J G xi (X)U p (X)dn(X), hf = J G x i(X)U h (X)dQ(X), 
hf(x) + hf(x)f(x), hf = j G xixj (X)U p (X)dQ(X), hf = J ' G xix i(X)U h (X)dQ(X), 



where all the /i's are known functions of x, and i, j = 1,2, Using this fact, the following 

types of relation (37) open different scenarios: 

1. F is an algebraic expression of its arguments, leading to an algebraic equation for /. The 
simplest case is, of course, that of a linear equation, leading to a linear algebraic equation 
for f(x). The simplest example of linear algebraic relation is obtained imposing that one 
of the fields be a "given function of the coordinates, interpretable as an external forcing" . 

2. F is a multidimensional linear partial differential equation of any dimension and order, 
either with constant or variable coefficients. This leads to a linear PDE for f(x) having 
the same dimension and order, and always variable coefficients (due to the functions h in 
equations (38)). 

3. F is a nonlinear PDE, whose dimension m is lower than the dimensionality n of the system 
of PDEs. This leads to a nonlinear PDE for f(x) in m dimensions, but with variable 
coefficients. In this case, the dressing procedure allows one to replace the nonlinear PDEs 
under investigation by a nonlinear PDE for f(x) in lower dimensions (a "reduction of 
complications"). 

Among all these cases, the most remarkable one is when F is a linear algebraic relation, 
since f(x) can be found explicitly and the solution manifold may be constructed analytically 
(see Sec. 4. 3 for more details on this point). 

The equation defining G(X;x), together with the relation (37) among the fields, provides 
the completeness of the nonlinear system of PDEs for the fields v, v 1 and w tj generated by the 
dressing. 

3.2 Examples 

In this section we consider some basic examples of partially integrable PDEs, corresponding to 
special definitions of G(X; x), and to particular relations (37). 

3.2.1 The simplest nonlinear partially integrable PDEs 

The simplest possible case corresponds to a function G independent of x. Then we have the 
only additional field w 00 , since w !J = 0, i,j > 0. We consider two examples of relation (37). 
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1. The relation (37) is chosen as follows: 



F : w (x) = exp 



E 

,i=i 



CLjX 



(39) 



where a,j are constant diagonal matrices. Then equation (36) with n = yields: 
A j = (a 2 - a ± B 2 + [B 2 , v])' 1 (aj - ai Bj + [Bj,v]). 



(40) 



Consequently, equations (33), for each particular choice of j, involve just the two fields v and 
v 1 . Thus we need two equations of this type to close the system, say j = 3,4, obtaining the 
following system of 2 matrix equations in 4 dimensions for the matrix fields v, v 1 : 



, -l 



(L 2 iv - [B 2 , (a 2 - a x B 2 + [B 2 , v}) 
(L 31 v - [B 3 , ui]) (o 3 - aiS 3 + [B 3 , v]) = (L 41 v - [B 4 , v^) (a 4 - a ± B 4 + [B 4 , v]) . 

2. The relation (37) is chosen as follows: 

F : v 1 (x) = 7(x), 

where j(x) is an arbitrary matrix function of x. In this case we choose: 

A\x) = {L 21 w 00 (x)y 1 L jl w 00 (x), j =3,4, 



(41) 



(42) 



(43) 



and we obtain the following system of 2 matrix equations in 4 dimensions for the two matrix 
fields v, w 00 : 



(L 21 v(x) - l2 (x)) (L^w 00 ^))- 1 = (L 31 v(x) - 73(a)) (L 3 W))" 1 
(L 4 it;(a:)-74^))(^4i^ 00 ^)) _1 , 

7j (x) = [B jn (x)l j = 2,3,4, 



(44) 



depending on the arbitrary forcing "f(x). 

In section Sec. 4 we will see that the space of analytic solutions of these two systems is 
2-dimensional. 



3.2.2 Partially integrable n-dimensional PDEs 

In this section we first derive nonlinear PDEs in 4 dimensions, and then we turn to nonlinear 
PDEs in arbitrary dimensions. 

Let G be defined by the following equations 

G xj = ajG x i, j = 1, . . . ,7i, (45) 

where aj are constant diagonal matrices. 

The special form (45) of G implies that the fields w j0 (x), w jn (x) can be expressed in terms 
of the fields w w (x), w 11 (x) respectively: 

wi°(x) = ajio 10 (x), w^ k {x) = ajakW 11 (x) j, k — 1, . . . , n. (46) 
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Then all the equations (33) and (36) reduce to the 6 equations (for j = 3,4): 



Ljiv - [B v v l ] = (L 21 v - [B 2 y\)Ai, (47) 
L,i«; 00 - a jW w + w^Bj = (L 21 w 00 - a 2 w 10 + w 10 B 2 )A>, (48) 
L jlW 10 - ajw 11 + w n Bj = (L 2l w 10 - a 2 w u + w n B 2 )A j , (49) 

Equations (49) can be viewed, for instance, as defining the matrix fields A j , j = 3,4; sub- 
stituting these definitions of A^ in equations (47) and (48), we obtain 4 matrix equations for 
the 5 matrix fields v, v 1 , w 00 , w 10 , w 11 . Before performing such operations, it is convenient to 
introduce a more compact notation defining the matrices 

E] := Lji(v) - [Bj, v 1 }, Ef° := Lj^w™) - a jW lk + w lk B j: j = 3, 4, k = 0, 1. (50) 

Then equation (49) yields 

A> = (Ef°)- 1 Ef°, j = 3,4 (51) 
and the nonlinear system (47,48) takes the form 



E^Ef)- 1 = El{Ef°)-\ j = 3,4, (52) 
EfiEf )' 1 = Ef\Ef°Y\ j=3,4. (53) 



This 4 dimensional system is not closed, consisting of 4 equations for the 5 functions 
v, v 1 , w 00 , w 10 , w 11 . This indeterminacy is consistent with the fact that the above equations 
are generated by the linear integral equation (6), which possesses a 1 dimensional space of 
homogeneous solutions. Therefore all the solutions of equations (52-53), constructed by the 
dressing procedure, contain an arbitrary function f(x). 

To fix this arbitrary function, we use the relation (37) which, for this example, reads: 

F(v , v 1 , w 00 , w 10 , w 11 ) — 0. (54) 

As we discussed above, this relation is largely arbitrary; the only requirement is that it should 
give rise to a solvable equation for the function f(x) (see Sec. 4. 3 for more details on this point). 

We remark that the arbitrary relation (37) can be chosen in order to put equations (52,53) 
into a differential polynomial form. If we choose, for instance, the bilinear relation 



„10 



F : Ef" = L 21 (w 10 ) - a 2 w n + w n B 2 = T, (55) 



where T is a constant matrix, multiplying equations (52,53) from right by EJ W , one transforms 
them into a differential polynomial form. In addition, if B 2 = 0, the bilinear relation becomes 
linear: w 1 ^ — a 2 w n = T. 

We have established that, due to the novel hypothesis dimker^ = 1, the dressing algorithm 
allows one to construct a system of partially integrable PDEs in multi-dimensions which includes 
one largely arbitrary relation among the fields. This is a novel and surprising feature of the 
theory. 

It is interesting to remark that equations (52-53) admit the reduction 

v = w 00 ; v i = w io j ( 56 ) 
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which follows imposing that G(X; x) = C(A; s), as a consequence of the choice otj = Bj. In this 
case, equations (52-53) reduce to the following two nonlinear PDEs 

E f°( E f°yi = Ef^Ef )' 1 = Ef^Ef )' 1 (57) 

for the three matrix fields w 00 , w w , w 11 , supplemented by the (largely) arbitrary relation 

F(w 00 , w 10 , w 11 ) = 0. (58) 

A simple example of linear relation (58) is, for instance, 

F : w n (x) = j(x), (59) 

where 7 is an arbitrary function; then the nonlinear system is given by the two equations (57) 
for the two fields w 00 ,w 10 , with the arbitrary forcing 7(2;) (it is the system (1-4) presented in 
the introduction with a different notation). 

As we shall see in section 4.4, the above nonlinear PDEs in 4 dimensions possess an analytic 
solution space of dimension 2. 

One can generalize the above construction, to generate partially integrable PDEs in n 
dimensions (n > 4), whose manifold of analytic solutions has dimension (n — 2). This higher 
dimensional generalization is associated with a more general equation for G. 

For instance, in order to derive partially integrable PDEs in 5 dimensions, one chooses 
G(X; x) to be defined by the following equations: 

2 

G x j = ^2a jk G x k, j > 2, (60) 
fc=i 

where are constant diagonal matrices. This special form implies that the fields w^ , w^ 2 , j > 
2 can be expressed in terms of the fields w t0 , w th , i,k — 1,2: 



w 



2 2 

vP k = Y,ot i8 w sk , j > 2. /,• 1.2. (61) 

s=l s=l 



Then the system (33,36) reduces to the following 12 equations: 

L jl v-[B j ,v 1 ] = (L 21 v-[B 2 ,v 1 ])A>, 3 = 3,4,5, (62) 
2 

L n w k0 - a H wki + wklB j = ( L 2iw k0 - w k2 + w kl B 2 )A j , j = 3, 4, 5, k = 0, 1, 2. (63) 



=1 



Following the previous procedure, we use equations (63) for k — 1 to define the matrices A\ 
and we substitute these definitions into equations (62) and (63) for k — 0. Defining the blocks: 



E) = L^-^V], j = 2,3,4,5, (64) 
2 

^ fe ° ee Lj^-^a^ + w 1 ^-, j = 3,4,5, k = 0,1, 2, 



i=i 

E^"" ee L 21 w™ -w k2 + w«'B 2 , k = 0,1, 2, 
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we obtain: 



A* = {E%»)- X E?°, J =3,4,5. (65) 

Substituting (65) into (62) and (63), we obtain the following closed system of 8 matrix PDEs for 
the 8 matrix fields v, Vi, w 00 , w w , w 20 , Wu, Wi 2 , w 22 , consisting of the following 5 basic equations 

E](EfY 1 = El{Ef°)-\ j = 3,4, (66) 
Ef(Ef)-i = Ef{Ef°)-\ j = 3,4,5, 

supplemented by two of the remaining equations (63), for instance, those for j — 3,4, k — 2: 

EfiEf )- 1 = Ef°(Ef°)-\ (67) 

and by one (largely arbitrary) relation between the 8 matrix fields: 

F(v, V!, w 00 , wio, w 2 o, w u ,w 12 , w 22 ) = 0, (68) 

which is introduced in the same spirit as before. Other equations of the system (62) and (63) 
can be treated as a symmetries of eqs. (66, 67). 

This procedure can be generalized to an arbitrary number n of dimensions, introducing the 
following equation for G: 

n— 3 

G xj = ^ a jkG x k, j > n - 3, (69) 
k=i 

where are constant diagonal matrices and n > 4. The resulting nonlinear PDEs, possessing 
the same block structure as their lower dimensional analogies, will have dimensionality n and, 
as we shall show in section 4.4, will be characterized by a manifold of analytic solutions of 
dimension n — 2. 

The possibility to increase the dimensionality of the PDEs and, at the same time, to increase 
proportionally the dimensionality of the manifold of solutions is due to the combined effect 
of the hypothesis: dim ker\I/ = 1 and of the introduction of the fields w tj . Indeed, i) the 
property (25) implies the multidimensional linear problems (30) and, via equations (51,65), the 
nontrivial mixing of the fields v, v 1 and w 1 ^; ii) the matrix fields are associated with the 
outer dressing function G, whose dimensionality can be increased without obstacles (see the 
(n — 3)-dimensional equation (69) ), while the matrix fields v, v 1 are associated with the matrix 
function C, an ingredient of the classical dressing method, whose dimensionality is severely 
constrained (see the 1-dimensional equations (10)). 

Note that, by construction, the derived systems possess higher symmetries: eqs. (52,53) with 
j > 5 and eqs. (66,67) with j > 6. 

We end this section elaborating on the dimensionality of the space of analytic solutions 
constructed in this section. Consider, as an illustrative example, the system of equations 
(57), (59) in n = 4 dimensions for the two matrix fields w 00 ,w w , and interpret x 4 as time 
variable. One can view the first equation E™ {E 2 ) _1 = Ef (Ef ) _1 as defining, in principle, 
w 10 in terms of w 00 and its partial derivatives w ^, j = 1, 2, 3. Replacing this relation into the 
second equation E™ (Eg y 1 = Ef (Eg )~\ one obtains a single equation for w 00 , depending 
linearly on w ® and nonlinearly on the other derivatives. Since, as we shall see in Sec. 4. 4, our 
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dressing algorithm generates analytic solutions depending on an arbitrary matrix function of 
n — 2 = 2 variables, the space of analytic solutions of (57), (59) has dimension n — 2 = 2. 
Similar arguments can be used for the systems (41), (44) and (66)-(68), which exhibit a single 
equation with first order derivative with respect to t — x n . Slightly different is a system (52)- 
(54), in which two equations involve first order ^-derivatives of the two functions v and w 00 . 
Accordingly, its solution space depends on two arbitrary functions of two variables, see (114) 
of 4.4. Thus we have established that the space of analytic solutions of all the examples of this 
section is (n — 2) -dimensional. 



3.3 Compatibility of linear spectral problems versus nonlinear PDEs 

It is well-known that integrable PDEs arise as the compatibility of overdetermined systems of 
linear problems for some eigenfunction U(X;x). For instance, the N-wave system (18) in 2+1 
dimensions is the integrability condition for the Lax pair (15). 

Such a picture is lost in our case, since there is no direct algebraic way to construct the 
nonlinear equations (32) as the compatibility condition of the linear systems (30). Indeed, the 
compatibility between equations (30) for j ^ k leads to equation 



U([B k ,L n v - [B j:Vl }} - [Bj,L kl v - [B k , Vl ]]) + (L 2 \U) (L kl Aj - L n A k )- 
(L 21 U) x i (A k Bj - AjB k ) + {L 21 U) xk A j - (L 21 U) x3 A k = 0, 



(70) 



from which one cannot infer anything, since the terms (L 2 iU) x j and {L 2 \U) x k are not inde- 
pendent matrix functions of A, being expressible in terms of L 21 U and U . But, to obtain such 
expressions, one has to use additional structure, i.e., the dressing equation ^Lj-JJ = 0. Differ- 
entiating it with respect to x k and using equation (8), one obtains the homogeneous equation 

* [(LjiU) x k + UB k (L jlV - [B J ,v 1 ])] = 0, (71) 

which, due to (25), implies that 

(L n U) xk = -UB k (L n v - [B vVl ]) + {L 21 U)P\ j + k, (72) 

where f nk are functions of x only. Substituting these relations (with j = 2) in (70), one obtains: 



U([B k , E] - E%Aj] - [Bj, El - E v 2 A k ]) + 

(Lai CO [L kl Aj - L n A k + f 21 (A k Bj - A 3 B k ) + f 2k Aj - f 2j A k ] = 0, 

and the independence of U and L 2 iU implies 

[B k , Ej - E^Aj] = [Bj, El - E%A k ], 

L k iAj — L jl A k + f 2 i(A k Bj — AjB k ) + f 2k Aj — f 2j A k = 0. 



(73) 



(74) 



We observe that equation (74a) does not imply directly the wanted equations (32). To obtain 
them, one should consider, instead, the compatibility between the linear problem (30) and 
equation (72), which leads to the following equation: 

UB k (E] - Elti) + (L 21 U) {A xk + f 2k A> - f k ) = 0. (75) 

Again the independence of U and L 2 iL implies the relations f jk = f 2k A j + A 3 k , together with 
the wanted equations (32). 

Summarizing, the linear 3-dimensional problems (30) contain only partial informations and 
their algebraic compatibility does not imply, alone, the nonlinear equations (32). 
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4 Solution space 



In the previous sections we have constructed partially integrable PDEs under the basic assump- 
tion that the integral equation (6) admits nontrivial homogeneous solutions. In this section 
we show i) how to choose the inner dressing functions in order to satisfy this assumption and, 
consequently, ii) how to construct the corresponding manifold of particular solutions of the 
partially integrable PDEs, expressed in terms of the dressing data $, C and G. Note that, 
in our case, the novel dressing function G appears, in comparison with the classical algorithm. 
The solutions of equations (9) and (10) are: 

$(A;x) = J %(X,k)e kB - x dk, (76) 

C(n;x) = J e« B - x C (q,fi)dq, (77) 

where k, q are scalar parameters and A, fi are vector parameters of dimension n — 3. Thus 
equations (8) yield: 



*(A,/x;*) = [ %(X,k)e^ B - x C (q,fi)^- + J:(X,f,), B X = I. 
J k + q 



(78) 



We remark that, in the case of .S-integrable equations, the integration constant S(A, fx) is chosen 
to be 8(X — /i). In our case, we need a special form for £ (see (82) below). 

It is quite standard to assume that the measure dQ(X) have support on an open domain 
V of the A-space, and on a disjoint discrete set of points D = {b±, . . . ,&m}, D n V = 0. 
Correspondingly, we use the following notation for the dressing functions. 

f 0(A;x)= / (f> {X,k)e kB - x dk, XeV, 
$(A;x) = <^ / (79) 

[ </> n (x) = U n0 (k)e kB - x dk, n=l,...,M, XeD, 

f c(X;x) = Je^ x c (q,X)dq, XeV, 
[ , ) \ c n (x) = fei B - x c n0 (q)dq, XeD, n—l,...,M, 



G(A;x) = { * G ^' 1 A/r (80) 

] g n {x), XeD, n = 1, . . . , M K ' 



m*) = { ^ n 1 »„ (81) 

1 u n (x), XeD, n = 1, . . . , M, v 1 



and we choose S(A, jj) in the form: 



5(\-fi), X,fieV, 
XeV n< 
XeD, n 

X,jieD, n, m = 1, . . . , M. 



S(A, fx) = { \ eV n ^ D ^ n= \- M M (82) 

' <r n (n), XeD, jJ,eV, n = l,...,M y ' 
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V 

M 



Then equation (6) reduces to the following system of M + 1 equations 

<t>{X-x) = [ M\k)e^ B -*c j0 (q)^ Uj (x)+ (83) 

O (A, k)e^ B - x c (q, fj,)u(fx; i)^^ + ]T a jW u j( x ) + AgD, 

+ q 3=1 

Mx) = £ / Mk)e^ B -*c j0 (q)^u j (x)+ (84) 
j= i J K + « 

| ^(Jb) e (^c b ( 9)M )«(^x)^M + 

M M „ 

y~](7 n jttj(x) + ^ / a n (fi)u(fi;x)d£l(/j,), n = l,...,M, 



for the unknown matrix functions u(X] x), X E V and Uj(x), j — 1, . . . , M. 

Once the solution is obtained, one constructs the matrix fields v, v 1 , w 1 ^ using equations 
(13, 17, 34): 



v(x 
v 1 (x 
w°°(x 
w j0 {x 
w lj (x 



M 



V 



c(A; x)u(X] x)dfl(X) + £^ Ck(x)uk(x), 

k=i 

M 

I [d xl c(X; x)]u(X; x)dft(X) + £ [d xi c k (x)]u k (x) , 

T"4 fe=l 

M 

#(A; x)u(X; x)dQ(X) + £ g k (x)u k (x), 

k=l 

/M 
[d Xj g(X; x)]u(X; x)dQ(X) + £ [d Xj g k (x)]u k (x), 
k=i 

M 

[d Xi d Xj g(X;x)]u(X;x)dtt(X) + £ [^^^(x)]^^) 



(85) 



77 



V 



k=l 



4.1 The condition dim ker\I/ = 1. 

Now we have to provide the condition drinker^ = I. We base our considerations on well known 
facts of the theory of linear integral operators. If the homogeneous equation 



(86) 



#(A, /I] x)H(/i] x)dQ(fj) = 
has a nontrivial solution, then its adjoint equation 

J H(\;x)V(\,ii; x)dQ(X) = 

16 



(87) 



has a nontrivial solution as well. If dimker^ = 1, then the solution spaces of both equations 
(86) and (87) are one dimensional. 
In our case, equation (87) reads 



x 



dx^^X-x^Cin-x') +E(A,//) 



dil(X), x' = x\ x ^ x v. (88) 



In view of the independence of $ and E, this equation is splitted into two equations: 

J H(X;x)$(X;x')dQ(X) = J H(X; x)E(A, /j)dQ(X) — 0, (89) 

which have to be satisfied for all x, x' and /i. This means that H is independent of x and, due 
to equation (76), the following two conditions must be satisfied: 

J H(X)%(X,k)dQ(X) = 0, Wk, (90) 

J H(X)E(X,n)dQ(X) = 0, V/i (91) 

for the existence of a nontrivial solution of the homogeneous equation (86). 

It is important to remark that, at the same time, the condition (90) provides also the 
solvability of the inhomogeneous integral equation (6). Therefore no further constraint must 
be imposed. 

We consider a particular way to satisfy conditions (90) and (91), choosing 

where the matrices Aj are constant and nonsingular, so that the conditions (90,91) are con- 
straints only for the discrete parts of $ and E: 

MM M 

A ih = A i°i M = = °' ^ G 13 ' n = l,...,M. (93) 

j=i j=i j=i 

Due to (93), we have only (M — 1) independent equations in the system (84) and, conse- 
quently, the solutions Uj(x) are constructed up to an arbitrary function f(x). 

We remark that terms containing aj, Oij and dj may disregarded in equations (83,84). 
Indeed, the terms containing Oj in equation (83) can be incorporated in the first term of the 
RHS. Similarly, the terms with Oij and <jj in the equation (84) can be incorporated in the first 
and second terms of this equation. Thus we set o,j = bj = = without loss of generality. 

4.2 Degenerate kernel. 

The system of linear equations (83,84), supplemented by the conditions (93), has a rich manifold 
of solutions. To construct explicit solutions, we choose, as usual, a degenerate kernel: 

M 

c (g, fi) = M9)c 2 j(^)- ( 94 ) 
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In this case, equations (83,84) reduce to the following linear system of M + M equations 

M M 

<t>n(x) = u nj (x)uj(x) + j> nj (x)uj(x) + u n (x), n = 1, ... ,M, (95) 

3=1 3=1 
M M 

4>n{x) = ^ p nj (x)Uj(x) + ^ pnj^Ujjx), 71 = 1, . . . , M, (96) 

3=1 3=1 

for the matrix fields Uj(x), Uk(x), j = 1, . . . , M, k = 1, . . . , M, where: 

u>k(x) = J c 2 fe(A)u(A; x)dQ(X), (97) 

and where the given coefficients v n j , v n j , p n j , p nj - , (fi n are defined in terms of the spectral func- 
tions: 

Pn,(x) = /0„o(^)e^ B -c jO (g)f^, /^(s) = /^^(fljg, (98) 

0n(s) = / c 2re (A)0(A;x)(ifi(A), nO (fc) = / c 2 „(A)</> (A, k)dQ(X). 

v 

This algebraic system is obtained, as usual, applying the operator J C2 n {p.;x)(Kt{p)- to (83). 

r> 

Having constructed, from (95), (96), the u s (x) and the u s (x), one obtains the eigenfunction 
u(X; x) via the formula: 

/M M 
O (A, k)e kB ' x dk - ^ p^(A; x)uj(x) - ^ /5,-(A; x)mj(x), (99) 

3=1 3=1 

where: 

P,(A; x) = J O (A, fc)e( fc+ ^c j0 (?)^, (100) 
Pi(A; x) = y O (A, fc) e (^>*^.( ? )^. (101) 

At last, one constructs the matrix fields v, v 1 , to*- 7 , solutions of the nonlinear PDEs of Sec. 3.2, 
from equations (85). 

4.3 Fixing the arbitrary function f(x) 

Due to the constraint (93), the solutions Uj and Uj of the algebraic system (95) depend linearly 
on an arbitrary matrix function f(x). Then, via (85), also the fields v, v 1 , depend linearly on 
this arbitrary function. Such an arbitrary function is completely fixed by the largely arbitrary 
relation among the fields. 
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(102) 



To be more concrete, let us illustrate all these facts in the simplest case: c (g, p) — 0, M — 2. 
The constraint (93) implies that p 2 j = —A 2 x A\p\j, j = 1,2; then, from the homogeneous 
version of (96): 

= p 11 (x)H 1 (x) + p 12 (x)H 2 (x), 
= p 21 (x)H l (x) + p 22 (x)H 2 (x), 

one verifies that the second equation is consequence of the first, while the first equation admits 
the solution 

H^x) = Pu 1 (x)f(x), H 2 (x) = - P T 2 (x)f(x) (103) 

depending linearly on the arbitrary matrix function f(x) (compare with (26)). The general 
solution of the inhomogeneous algebraic system (96) is then given by: 



iii(a;) = 5Pn(aO0i(aO + Pu(x)f(x), 
u 2 (x) = Ipw( x )M x ) ~ Pii( x )f( x ), 



(104) 



depending on the arbitrary matrix function f(x) in a linear way as well (compare with (27)). 
Consequently, such a linear dependence on f(x) will appear, via (85), also in the matrix fields 



v, v 1 , w l K 



We remark that one could always identify f(x) with one of the u/s, say, with ui(x), obtaining 

ui(x) = f(x), u 2 {x) = Pi 2 (x) [<t>i(x) - p u (x)f(x)] ; (105) 

this identification, which clearly leads to a less symmetric formula than (104), seems to become 
more convenient when M > 2. 

Now we show, always in the simplest case: co(q, p) = 0, M = 2, how the arbitrary function 
f(x) gets fixed imposing the relation (37), which we choose in one of the forms pointed after 
equation (38). 

Equations of Sec. 3. 2.1. For the equations of this section, choose g(X) = 8(X — a). 
Then, using equation (85) for w 00 , the constraint (39) becomes the linear equation 



ho(x) + [h l (x)p l ^(x) + h 2 {x)p 12 (x))] 



4>i jx) 



(106) 



[h( x )pn(x) + h 2 (x)p 12 l {x))] f(x) = exp (£" =1 atf) 
for f(x), where 

hj{x) = - J g{\; x)p 3 {\- x)dQ(X) + 9j (x), j = 1, 2, h (x) = J g(X; x)0(A; x)dQ(X) (107) 



v 



V 



yielding the following explicit formula for /: 

f(x) = [h!(x)p^(x) + h 2 (x)p^(x))] 1 
h (x) - [hi(x)p^(x) + h 2 (x)p^ 2 l (x))] 



ex p (Er=i a i xi ) - 



(108) 



Analogously, using equation (85) for v 1 , the constraint (42) becomes a linear equation for f(x), 
whose explicit solution is 



f(x)=[h 1 (x)p u 1 (x) + h 2 (x) Pu 1 (x))]- 1 



li x ) - [hi{x)pn{x) + h 2 {x)p 12 \x))] 



(f>i(x) 



, (109) 
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where now: 

hj(x)= j qe* B - x c jQ (q)dq, J = 1,2. (110) 
Equations of Sec. 3. 2. 2. In this case, equation (69) implies that 



G(X;x) = J 



cxp 

V 




G {X',X)dtt{X'). (Ill) 



If, in particular, n = 4, as in equation (45), then the reduction G = C is admissible, 
identifying ctj = aiji = Bj and G = C (see (77)), and one obtains the nonlinear PDEs (57), 
supplemented by the relation (68). 

If this relation is an arbitrary linear relation between fields, it leads to a linear equation 
for f(x). For instance, if we choose (59), then, using the equation (85) for it; 11 , the constraint 
becomes a linear equation for f(x), whose explicit solution is 



f{x) = [h^p^ix) + h 2 (x)p 1 i(x))] [t(x)- 
Jg xlxl (X;x) ( f>(X;x)dQ(X)-[h 1 (x)p^(x) + h 2 (x)pT 2 1 (x))} *p 



(112) 



v 

where: 

hj(x) = - J g x i x i(\;x)p j (\;x)dtt(\) + g jxlxl (x), j = 1,2. (113) 
v 

In a similar way, one can treat more general relations in higher dimensions n. 
4.4 Dimensionality of the solution space 

The dimensionality of the space of analytic solutions generated by our dressing scheme is 
essentially defined by the dimensionality of two expressions 

fC(\;x)Q(\;x)d£l(\), j G(X;x)$(X;x)dtt(X). (114) 

The first term, involving the "inner" dressing functions C(X;x) and <3>(A;x), appears also in 
the classical dressing; its dependence on the space-time coordinates is severely constrained and, 
consequently, it carries dimensionality 2 (see (24) and the considerations made there). The 
second term involves the "outer" dressing function G(X;x), a novel feature of our dressing 
procedure; its dependence on the space-time coordinates is instead largely arbitrary, playing a 
crucial role in increasing the dimensionality of the solution space through the following novel 
mechanism. 

As we have seen in Sec. 4. 3, from the largely arbitrary relation among the fields one construct 
f(x) in terms of the spectral representations of the fields v,v l ,w lj . If such relation involves 
the fields w l \ whose spectral representations involve expressions like (114b), then the dimen- 
sionality of f(x) is not severely constrained. Since the matrix fields i>, v 1 , u> u , solutions of our 
PDEs, depend linearly on f(x), their analytic solution space is not severely constrained too. 
Using this argument, it is possible to establish easily the dimensionality of the space of analytic 
solutions of our PDEs. 

In Sec. 3. 2.1, the first term has higher dimensionality than the second, since G does not 
depend on x; it follows that the dimensionality of the space of analytic solutions of the 4 
dimensional PDEs constructed there is 2, like for integrable PDEs in 2+1 dimensions. 
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In Sec. 3. 2. 2, the outer dressing function is (111) and the second term in equation (114) 
reads (choosing G (X', A) = 8(X' — A)): 



cxp 



v 



5> K+ E 



k=n-2 



$ (A, k) exp [fc (£ • x)]dkdQ(X). 



(115) 



Since A is a vector parameter of dimension n — 3, and A; is a scalar parameter, the above 
expression has dimension n — 2, being parameterized by the arbitrary function $o(A, fc) of 72,-2 
variables. Consequently, n — 2 is the dimension of the constructed /(#), and of all the fields 
appearing in the nonlinear PDEs. 

Summarizing, the solutions we constructed depend on an arbitrary function of n— 2 variables 
(114b) and on an arbitrary function of 2 variables (114a). Then, in the exceptional case n — 4, 
the solutions depend on 2 arbitrary functions of 2 variables. This conclusion is valid for all the 
examples presented in the paper. 



5 Conclusions 

In this paper we have generalized the dressing method to construct systems of nonlinear PDEs 
in n dimensions (n > 3) i) possessing a manifold of analytic solutions of dimension n — 2 (a very 
large, but not complete, manifold), and ii) possessing higher symmetries. But the constructed 
PDEs do not seem to be the compatibility condition for overdetermined systems of linear PDEs, 
a characterizing feature of completely integrable systems in lower dimensions. 

The above properties indicate that they are examples of partially integrable PDEs in multi- 
dimension possessing a very large, but not complete, space of analytic solutions. 

A natural generalization of the algorithm presented in this paper consists in studying the 
case in which the integral operator ^ of the dressing problem exhibits a higher dimensional 
kernel: 

dimker* = D ker > 1. (116) 
In this case, equation (30) is replaced by 

D ker +1 

L ml U(X-x)+ (L nl U(X;x))A mn (x) = 0, m > D ker + 1 , (117) 

71=2 

and one needs D ker conditions on U to define the functions A mn . The study of the structure 
of the associated partially integrable equations, and of the dimensionality of the associated 
analytic solution space is postponed to future investigations. 
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